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EFFECTS OF CORRELATION ON THE STRUCTURE OF IMPURITY BANDS 

by Gust Bambakidis 
Lewis Research Center 

SUMMARY 

The effects of electron-electron interaction on the electronic density of states has 
been considered for a dilute system of scattering centers distributed randomly in an 
electron gas of moderate density. Using as the starting point an independent-electron 
t-matrix method, a modified t-matrix in the presence of the electron-electron inter- 
action is obtained. To illustrate the effects of correlation, the formalism has been 
applied to the simple case of a linear chain of attractive 5 -function potentials distributed 
randomly. The electrons are assumed to interact by means of a repulsive 5-function 
potential. The density of states in the impurity band is evaluated and compared with the 
results for independent electrons. It is found that the inclusion of electron -electron 
interaction produces a significant broadening and downward energy shift of the impurity 
band. 


INTRODUCTION 

Many papers dealing with the electronic structure of disordered systems have 
appeared. Because of the difficulty of the subject, for which there exists no general 
one-electron theory, these treatments have, for the most part, considered the electrons 
as noninteracting. Wolff (ref. 1) has treated the modification of a free-electron band in 
the presence of both the electron -impurity and the electron-electron interactions. How- 
ever he considers the case of very high impurity and electron densities (relative to the 
effective Bohr radius of the impurity), so that no localized states arising from single- 
impurity bound states can occur. The transport properties of interacting electrons in a 
system of random impurities have also been treated extensively by Langer (refs. 2 to 5), 
by Betbeder-Matibet and Nozieres (ref. 6), and more recently by Sigel and Argyres 
(ref. 7) and by Sigel (ref. 8). Here we describe a calculation of the effect of particle- 
particle interaction on the density of single -particle states in the band which arises from 
a bound state of an isolated impurity. In order to describe the impurity band without 
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many -particle effects, the t-matrix approach of des Cloizeaux (ref. 9) is used. This 
approach is based on infinite -order perturbation theory, and, consequently, many -body 
effects can be included in a straightforward way using standard propagator formalism 
(refs. 10 and 11). 


CALCULATION 

The Hamiltonian for a system of fixed-volume n at T « 0 containing N g particles 
interacting with N a impurities distributed at random is written. The impurities are 
assumed to be infinitely massive. Then 

Ng N e N e Ng Ng 

i=l i,j=l 

M3 

where 
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i=l 


Y/ V oa ' 2 Y V oe ^ 

i=l * i=l 



is the kinetic-energy operator for the i particle and 
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V a<£i> = Y V a ( £i " 
a - 1 


( 2 ) 


( 3 ) 


th 

is the potential energy of the i particle due to the interactions v (r. - R ) between the 
i particle and the impurities distributed among the sites R . All physical quantities 
are to be averaged over the (random) distribution of the R . We consider the system as 
being "neutral; ” that is, assume that there are potential sources which compensate for 
the potential energy of a particle in the average potential fields of the impurities V 

Ocl 

and of the particle gas V Qe . The perturbation calculation is carried out by using the 
second -quantized form of H: 
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The effect of the constant term -V Qe in equation (1) is to remove the diagonal elements 
<k| v e |k> from the particle -particle portion of H. Hence, k ^ k’ in equation (1’)- The 
constant term -V_„ in equation (1) simply shifts the scale of single -particle energies: 

Ud 

E = E’ + V oa = E , + N a < k|v a |k > 

where E 1 is the energy referred to the unshifted energy scale. (However, when v„ is 

a. 

long-ranged (e. g. , Coulombic), this energy shift is divergent and -V „ has to be in- 
eluded implicitly by removing the diagonal elements from the particle -impurity inter- 
action. ) 

We wish to obtain the single -particle propagator G(k, co) for the system. The density 
of single -particle states will then be given by 


D(E) = I Im ^ G (k, E) 
71 k 


E < E-, 


- Im V G(k, E) 
77 k 


E > E, 


(4) 
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(Unless otherwise indicated, the index k includes implicitly the spin quantum number 
a = ± 1/2.) First, consider the case without particle -particle interaction (v = 0). The 
Hamiltonian becomes a sum of N independent one-particle Hamiltonians H. = T. + 

V a (£j) > and it is sufficient to obtain G(k, oo) for the system of one particle plus impurities. 
The propagator can be written as a summation of an infinite series of multiple-scattering 

.r 

diagrams involving the free -particle propagator G Q (k, oo). Denoting G(k, oo) by k, oo 
and G Q (k, a>) by k , oo | yields 
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Next, the potential v„ is replaced by the particle -impurity t-matrix t (co), which is 

ci Si 

the sum of all multiple scatterings of the particle by the impurity a: 
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In reference 9, des Cloizeaux has introduced the notion of articulation points in these 
diagrams. By definition, an articulation point is associated with an interaction of the 
particle with the impurity ot if all scattering centers encountered by the particle before 
this interaction 0(j3 ^ a) are different from all scattering centers encountered after this 
interaction y(y ^ a). As a consequence of the averaging process, each diagram can then 
be reduced to a product of irreducible parts by cutting the propagator lines connected to 
each articulation point in the diagram. Propagation between articulation points can be 
described by a "self -propagator” S(k, oj), with matrix element 


<k 2 |S(k, cu) |kj) 


k, cj 


k 
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Upon associating the operator P a t a (o;) s N a /fi t a (o)) with each articulation point denoted 
by O, the diagrammatic expansion of G(k, w) can be redrawn as 
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leading to 


G(k, u>) = G Q (k, o>) + G a (k, c o) < k|P a t a (co) [l - P a S(k, co) t a (w)J 1 |k> 


(5) 
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The self -propagator S can, in turn, be expanded in powers of the impurity density with 

2 

the leading term, of order p„, containing the contribution of scattering by pairs of im- 

a 

purities. 

In considering the effects of electron-electron interaction, we are particularly in- 
terested in the broadening and shift of the impurity band at moderate to high electron 
densities. Consequently, the effect on a particle of a single impurity is examined in the 
presence of all the other particles. This approach leads to a modification of the particle - 
impurity t-matrix while neglecting correlation effects on the scattering of the particle 
between different impurities. For each ’’skeleton” particle -particle scattering diagram 
contributing to the propagator in the absence of this impurity, the corresponding diagram 
in its presence is obtained by replacing everywhere the free -particle propagator with the 
’’dressed" propagator appropriate to the system of a particle plus one impurity: 


k, o o 
k' u> 




k, a) 



k\ co 


Consider irreducible diagrams: that is, those which cannot be cut into two others by 
cutting a dressed propagator line. We denote by S(k,kJ co) the sum of all irreducible 
diagrams, taking the particle from state k' to state k, with energy Hcu: 



Because the presence of the impurity disrupts the homogeneity of the particle gas, the 
term S(k,k' co) differs from the usual self-energy in that it is nondiagonal. For this 

reason, Dyson’s equation for the modified propagator | for the system of particle 
gas plus one impurity cannot be summed into closed form. We can write only 
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The effective t matrix is defined in terms of 

a 


by 


E, “ 

- k . 

k, w I 

where k, c o is the propagator G e (E, u>) in the presence of the particle gas only, 

eff 

To obtain t we must know 2. To make the situation tractable let us in fact approxi- 
mate 2 by 6^ j f t2 e (k, u>), where 2 g is the (diagonal) self-energy without the impurity. 

The expression for ft then becomes 
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where 




+ 


k, co 


k, co 



+ 



G Q (k, a >) 

1 - S e (k, w)G 0 (P, cu) 


( 8 ) 


Let us now consider the simple case of a disordered linear chain of attractive 6- 
function potentials. In this case, 

V a (x - x„) = -v*5(x - X a ) 


and the t -matrix is 


-v 


<k|t (w)|k* > = 


n 


iv. 


2(w + io) 


1/2 


(a < 0 


in atomic units (R = 1, m„ = 1/2, e /4 ttf = 2). This has a singularity at the bound-state 

a 9 e ® 

energy = -(v^/2) of a particle in the field of one impurity. The fact that the matrix 

element of t is independent of initial and final momenta simplifies matters considerably 
a 


since equation (7') for 


G a-e^ k,k ’’^ H 


k,cu A 


k, u> 


the propagator for the system of 


particles plus one impurity, can be summed into closed form: 
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G a _ e (k,k/co) = G e (k,a>) • 5^, + G g (k; w)t a (w)G e (k, w) 


+ G g (kj cu)t a (aj)^ Yj [ G e ( k l’ " G 0 ( k i> w )] 

ki 


■V“) G e^» w ) 


+ • • • 


G e ( k , w) • 6^ k , + G e (k;w)<j ^ t^ +1 (w)| n (w)j>G e (k, w) 

i n=0 


= G e (k, 10 ) • 6 k ^ k , + G e (k| w)t a (w) [l - t a (w)J(w)] 1 G e (k, w), 


where 


J(w) = Yj [ G e ( k i> w ) “ G 0 ( k i» w 3 

ki 


(9) 


G Q (k 1 , w)S e (k 1 , co)G 0 (k 1 , a>) 

1 - S e (k 1 ,a))G 0 (k 1 ,aj) 

and t (w) has been written for < kit (w) J k T >. From equation (6) we see that 

tf(co) = [l “ s e ( k ’ w)G 0 (k! co)] 1 t a (w) [l - t a (u>)S(o>)] _1 • [l - S e (k,co)G 0 (k, w)] , 

and is dependent upon k' and k. To facilitate comparison with the results using the 
bare t -matrix, this dependence is ignored: 

tfM « t a (w) [l - t a (w)J(w)] _1 . (11) 

The quantity £>{<x>) has been evaluated when the particles or "electrons” interact by 
means of a repulsive 5 -function: 



(10) 
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v e ( x i " x j) = v o 6(x i - x j)' 


( 12 ) 


The self-energy was calculated in the ladder approximation to lowest order in the 
particle -particle t-matrix t e (q, £): 


s e (k, oj) « r q - IT, £ - w j^) q - k, £ - 

te^’ ^ • x — > f°° r _ _ _ -i 

+ z rr.IZ r = - • 2 ^ J d£ |< k |t e (q, £) |k > - < k|v e |k>jG 0 (q-k,£-cu) 

q-k, £-w q 

d'£ < q - k |t e (q, £) |k > G Q (q-k, £-uj) (13) 

q 


In this case < k|t e (q, £) |k' > depends only on q and £ and can be obtained in closed 
form as 


t e (q, £) = 


1 - 


o lim 


n 
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£ - k 2 - (q - k) 2 + i5 
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£ - k 2 -(q - k) 2 - 15 


(14) 


DISCUSSION 


The expression which results when the free -particle propagator G Q (k, cu) appearing 
in equation (5) is replaced by G g (k, co) and the bare particle -impurity t-matrix t a by 


t 


eff 


is 


11 



G(k, w) = G e (k, w) + Gg(k, co) <k jP a tJ ff (cu) [l - P a S(k, w)tj ff (u)] _1 |k >. (15) 

This equation reduces to the correct limit when either v 0 or v vanishes. Upon com- 

cff ” ^ 

puting G from equation (8) and t from equation (11), the resulting density of states 
(per unit length) obtained for negative energies is shown in figure 1, together with the 
result for independent electrons. In both cases we have used v a = 2Ry - a and P„ = 

1 O i O cl 

1/3 a” , while in the former case we have taken v e = v a and P = 1 a” 1 . The lowest 
o o o e o 

order expansion in P of the self -propagator was used. (The self -propagator S(k, oj) 

a. — 

consists of the sum of two terms A(a>) + B(k, co). To lowest order in p„, only the first 
term has been retained as giving, in this particular case, the dominant contribution to 
D(E) in the neighborhood of the maximum (see ref. 9).) Also shown is a machine calcu- 
lation for independent electrons, reported in reference 6. We see that the effect of 
particle -particle interaction on the impurity band is to displace it to lower energies as 
well as to broaden it. For the parameter values used here the full width at half maximum 
with particle interaction (~0. 42 Ry) is comparable to that occurring in the machine cal- 



Energy, E, Ry 


Figure 1. - Single-particle density of states 
scatterers. Scatterer density, p a = 1/3 a~l 


per unit length, in the impurity band for a linear chain of random 
; particle density, p e = la" 1 . 
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culation for independent particles (~0. 30 Ry). To estimate the shift of the impurity band 
with respect to the Fermi level, the following approximate relation for a filled band 
(P e > 2P a ) was used in obtaining Ep,: 



[d(E) - D 0 (E)] dE ~2P £ 


(16) 


or 



I D (E)dE ~P 
Si e 


2 P, 


(16’) 


where D g (E) is the contribution to the state density from the first term in equation (15). 
Equation (16) is a good approximation in the independent -particle theory, for which 

(l/f2)D(E)dE ~ 2 p In general, we would expect it to be applicable in the tight- 

binding limit (P a small relative to the bound-state radius of an isolated impurity) for 

v small compared to v„. Using equation (16') we obtain E-,- ~ 0. 27 Ry, eL — 0.03 Ry, 
e a r r 

so that the shift is 


(Eg - Eg) - (E B - E p ) ~ -0. 49 Ry 


The integrated state density up to E = Eg is 



— D(E)dE = 0. 737 

n 


states ^ p 
a o 


Thus equation (16) underestimates the shifted Fermi level in this case (v = v„). 

The above results indicate that while broadening of the impurity and conduction bands 
due to interparticle interaction increases their overlap, there is a net downward shift 
of the center of the impurity band with respect to the Fermi level. 

In three dimensions, the effect of particle -particle interaction on an impurity band 
should be more pronounced than in one dimension. In three dimensions, strong screen- 
ing of an impurity at high electron densities can cause it to ionize completely, whereas 
a one -dimensional potential well always possesses at least one bound state. We believe 
that the overall results of this calculation are reasonable. In treating a more realistic 
(e. g. , three dimensional) case, however, the computational problem becomes rather 


13 


formidable. An accurate determination of 2 e (k, u>) for the Coulomb interaction for a 
wide range of c o, for example, is difficult at intermediate to low electron density, where 
correlations are important. Such a calculation would be of interest in any situation where 
one probes the state density arising from tightly bound outer electron states of impurities 
and/or atoms of a disordered host. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 28, 1970, 

129-03. 
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